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Abstract 

The well-known algorithm for summing divergent series is based on the Borel 
transformation in combination with the conformal mapping. A modification 
of this algorithm allows one to determine a strong coupling asymptotics of 
the sum of the series through the values of the expansion coefficients. An 
application of the algorithm to the /3-function of <£> 4 theory leads to the 
asymptotics 0(g) = (3oog a at g — > oo, where a ~ 1 for space dimensions 
d = 2,3,4. The natural hypothesis arises, that the asymptotic behavior 
is (3(g) ~ g for all d. Consideration of the "toy" zero-dimensional model 
confirms the hypothesis and reveals the origin of this result: it is related to 
a zero of a certain functional integral. A generalization of this mechanism 
to the arbitrary space dimensionality leads to the linear asymptotics of /3(g) 
for all d. The same idea can be applied to QED and gives the asymptotics 
/3(g) = g, where g is the running fine structure constant. A relation to the 
"zero charge" problem is discussed. 



1. Introduction 

It is commonly accepted that summing divergent series can give important and non- 
trivial information. It will be demonstrated below that sometimes we can obtain even more: 
summation of the series allows to guess the exact result and then this result can be proved. 

Our main interest is a reconstruction of the Gell-Mann - Low function (3(g) for actual 
field theories from its divergent perturbation expansion. We describe the summation pro- 
cedure in Sec. 2 and illustrate it for the case of ip 4 theory in Sec. 3. The arising hypothesis 
on the linear asymptotics (3(g) oc g is tested in Sec. 4 in the zero- dimensional limit, while 
Sec. 5 gives its justification for any dimension d < 4. The same idea is applied to QED in 
Sec. 6. Finally, Sec. 7 discusses some problems arising in relation to the obtained results. 



2. Summation procedure 
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Let us consider the typical problem in field theory applications. A certain quantity 
W{g) is defined by its formal perturbation expansion 



W(g) = J2 W N (-g) N (1) 

N=0 

in the powers of the coupling constant g. The coefficients Wn are given numerically and 
have the factorial asymptotics at N — > oo, 

W a N s = ca N T(N + b) , (2) 



which is a typical result obtained by the Lipatov method |15J . One can see that the 
convergence radius for (1) is zero. The problem arises, can we make any sense of the series 
(1) and find W(g) for arbitrary g. 

The conventional treatment of the series (1) is based on the Borel transformation 



W(g) = I dxe~ x x h °- l B{gx) , (3) 



B(z) = £ B N (-z) N , B N = —^— , (4) 



N=0 



T(N + b ) 



relating the function W(g) with its Borel transform B(z), while B(z) is given by a series 
with a factorially improved convergence; bo is an arbitrary parameter, which can be used for 
optimization of the procedure. Under the proper conditions, Eq.3 is an identity obtained by 
interchanging of summation and integration and using a definition of the gamma-function. 
In the general case, Eqs.3,4 give a definition of the Borel sum for a series (1). In what 
follows, we identify the function W(g) with the Borel sum of its perturbation series. In 
the case of ip A theory, it is possible to test a validity of such identification in one and zero 
dimensions [21] and to prove the Borel summability in two and three dimensions [To^ Ej. 

It is easy to show that the Borel transform B(z) has a singularity at the point z = — 1/a 
(Fig. l,a) determined by the parameter a in the Lipatov asymptotics (2). The series for 
B(z) is convergent in the disk \z\ < 1/a, while we should know it on the positive semi- axis, in 
order to perform integration in the Borel integral (3); so we need an analytical continuation 
of B(z). Such analytical continuation is easy if the coefficients Wn are defined by a simple 
formula, but it is a problem when they are given numerically. 

The elegant solution of this problem was given by Le Guillou and Zinn- Justin in 1977 
[8J. It is based on the hypothesis that in field theory applications all singularities of B(z) He 
on the negative semi-axis. This hypothesis can be proved in the case of <p* theory [211.111 If 



1 A validity of this hypothesis is frequently questioned in relation to possible existence of the renormalon 
singularities [7]. Such singularities can be easily obtained by summing some special sequences of diagrams, 
but their existence was never proved, if all diagrams are taken into account [2j. The present results for the 
asymptotics of the /3-function (Sees. 5, 6) are in agreement with a general criterion for absence of renormalon 
singularities [55] and a proof of their absence for ip 4 theory [31] (see a detailed discussion in [55]). 
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Figure 1: (a) The Borel transform B(z) is analytical in the complex plane with the cut 
(— oo, — 1/a); (b) Its domain of analyticity can be conformally mapped to a unit disk in the 
u plane; (c) If analytic continuation is restricted to the positive semi-axis, then a conformal 
mapping can be made to any domain, for which the point u — 1 is the nearest to the origin 
of all boundary points; (d) An extreme case of such domain is the u plane with the cut 
(l,oo). 



3 



such analytical properties are accepted, we can make a conformal transformation z = f(u), 
mapping the complex plane with the cut (Fig. l,a) into the unit disk \u\ < 1 (Fig. 1,6). If 
we re-expand B(z) in the powers of u, 

OO oo 

B(z) = £ B N (-z) N y = f{u) B(u) = £ U N u N , (5) 



N=0 



N=0 



then such series will be convergent for any z except the cut (— oo, — 1). Indeed, all singular 
points P, Q, R, . . . of B(z) lie on the cut, and their images P, Q, Q', R, R', . . . in the u 
plane appear on the circle \u\ = 1. The re-expanded series in (5) is convergent for u lying 
within the unit circle, but the interior of the circle \u\ < 1 is in one-to-one correspondence 
with the analyticity domain in the cutted z plane (Fig. 1, a). 

Such conformal mapping is unique (apart from trivial modifications), if we want to 
make an analytical continuation to the whole domain of analyticity. In fact, such strong 
demand is not necessary since we need B(z) only at the positive semi-axis, in order to 
produce integration in (3). If we accept that the image of z = is u = and the image 
of z — oo is u — 1, then we can make a conformal mapping to any domain, for which the 
point u — 1 is the nearest to the origin of all boundary points (Fig. l,c). The series in u 
converges for \u\ < 1, and in particularly at the interval < u < 1, which is the image of 
the positive semi-axis. 

The advantage of such conformal mapping consists in the possibility to express the large 
g asymptotics of W(g) in terms of the expansion coefficients Wn- Indeed, the divergency 
of the series in u is determined by the nearest singular point u — 1, which is an image of 
infinity: so the large iV behavior of the expansion coefficients U N is related to the strong 
coupling asymptotics of W{g). In order to diminish influence of other singular points 
P, Q, Q', . . ., it desirable to remove these points as far, as possible. Thereby, we come to 
an extremal form of such conformal mapping, when it is made on the whole complex plane 
with the cut (1, oo) (Fig. l,d). Mapping of the initial region (Fig. l,a) to the region of 
Fig. 1, d is given by a simple rational transformation 



u 



a(l — u) 

for which it is easy to find the relation of Un and B^, 



(6) 



N p> 

U = B , U N =J2 -f(-l) K C«ll (N>1), (7) 



K=l a 



where = N\/K\(N — K)\ are the binomial coefficients. If W(g) has a power law 
asymptotics 

W(g) = W oo9 a , g^oo, (8) 
then the large order behavior of Un 

U N = U^N a - 1 N -> oo , (9) 
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00 a^T(a)T(b + a) 1 j 

is determined by the parameters a and W^. Consequently, we come to a very simple 
algorithm [21]: the coefficients Wn of the initial series (1) define the coefficients Un of re- 
expanded series (5) according to Eqs. 4, 7, while the behavior of Un at large N (Eqs. 9, 10) 
is related to the strong coupling asymptotics (8) of W(g). 

If information on the initial series (1) is sufficient for establishing its strong coupling 
behavior (8), then summation at arbitrary g presents no problem. The coefficients Un are 
calculated by Eq.7 for not very large N, and then they are continued according to their 
asymptotics (9). Consequently, we know all coefficients of the convergent series (5) and it 
can be summed with the required accuracy. 

Few comments should be made to avoid a misunderstanding. The conformal mapping 
corresponding to Fig. 1, b provides (for fixed z) the fastest convergence rate for the u series 
[5J, and is cited as "optimal" in the literature. It may look preferable to use this algorithm 
and extract the asymptotics of W(g) from the summation results. In fact, all investigators 
of the strong coupling region [11] [131 EH1 El] independently came to the same conclusion 
that the asymptotics of W(g) should be estimated before any summation!! On the other 
hand, the fastest convergence is a distinctive excellence only if Wn are known exactly. In 
the presence of round-off errors, the uncertainty in Un grows as 5.8^ for Fig. 1,6 and as 
2 N for Fig. l,d [21]; more than that, the latter (but not the former) algoritm is stable in 
respect to smooth errors (like interpolation ones) |24j . and it has a crucial significance for 
the following applications. 



3. Application to tp A theory 

The described algorithm was successfully tested for a lot of simple examples [21], and 
now we can apply it to a reconstruction of the Gell-Mann - Low function (3(g) of quantum 
field theories. This function enters the Gell-Mann - Low equation which describes the 
behavior of the effective charge g as a function of the length scale L: 

(I<J M • (ii) 



d\nL 



The most interesting problem is an appearance of the /3-function in relativistic theories, 
like four- dimensional ip 4 theory or QED. In this case, the expansion of (3(g) begins with 
the positive quadratic term and the effective charge g grows at small distances|f| (Fig. 2): 
it is interesting to find the law of this growth in the strong coupling region. 

According to the classification by Bogolyubov and Shirkov [3], there are three qualita- 
tively different possibilities (Fig. 3): (1) if /3(g) has a zero at some point g*, then the effective 



2 For example, it is clear from the described algorithm, that one cannot find a correct asymptotics of 
W(g), if he does not know a correct asymptotics of Un- 

Equation (11) is valid for L<to _1 , where to is a mass of the particle; in the region L>to _1 , g remains 
constant and equal to its observed value g b s ■ 
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Figure 2: Effective coupling g as a function of the length scale L in four- dimensional y? 4 
theory and QED. 




Figure 3: Three qualitatively different situations according to the Bogolyubov and Shirkov 
classification. 
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coupling g tends to g* at small L; (2) if /3(g) is non-alternating and has the asymptotic 
behavior g a with a < 1, then grows to infinity; (3) if non- alternating 0(g) behaves 
at infinity as g a with a > 1, then g(L) is divergent at some finite L$ and the dependence 
g(L) is not defined at smaller distances: the theory is internally inconsistent and a finite 
interaction at large distances is impossible in the continual limit. To distinguish between 
these three possibilities, one needs to know the /3-function at arbitrary g, and in particular 
its asymptotic behavior for g — > oo. 

One can attempt to solve this problem by summation of the perturbation series, 

0(g) = f3 2 g 2 + (3 3 g 3 + ... + (3 L g L + ... + ca N T(N + b)g N + ..., (12) 

having in mind that several first coefficients (till (3l) are known from diagrammatic calcu- 
lations and their large order behavior is given by the Lipatov method. The intermediate 
coefficients can be found by interpolation, the natural way for which is as follows. It can 
be shown that corrections to the Lipatov asymptotics has a form of the regular expansion 
in 1/JV: 

N = ca»r{N + b){l + ^ + £ + ... + £L + ..) . (13) 

One can truncate this series and choose the retained coefficients Ak from correspondence 
with the first coefficients 02i ■ ■ ■ j 0L\ then the interpolation curve goes through the several 
known points and automatically reaches its asymptotics. To variate this procedure, one 
can re-expand the series (13) in the inverse powers of iV — N, 

(3 N = ca N T(N + b){l + — A * +...+ Ak +.\ , (14) 
V 1 \ N-N (N-N) 2 (N- N) K J V 1 

and obtain a set of interpolations, determined by the arbitrary parameter N. 

In the case of four-dimensional <£> 4 theory, a realization of this program [21] gives the 
non-alternating /3-function (Fig. 4, a), with the results for the exponent a shown in Fig. 4, b. 
The exponent a is practically independent on N, and only its uncertainty depends on this 
parameter. If we take the result with the minimal uncertainty, we have a value a = 
0.96 ± 0.01, surprisingly close to unity0 

Something close to unity is obtained also in two and three dimensions [TSl [TTJ] (Fig. 5). 
The natural hypothesis arises, that (3(g) has the linear asymptotics 

0(9) ~ 9 , 9 -> 00 (15) 

for arbitrary space dimension d. If this hypothesis is correct, then there is a natural strategy 
for its justification: 

(i) to test it in a simple case d = 0; 

(ii) to find out the mechanism leading to this asymptotics; 



4 Estimation of errors was made in a framework of a certain procedure worked out in [24] . Subsequent 
applications have shown that such estimation is not very reliable. 
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Figure 4: (a) General appearance of the /3-function in four-dimensional y? 4 theory according 
to [21] (solid curve), and results obtained by other authors (upper, middle, and lower dashed 
curves correspond to [HI [131 EH] respectively), (b) Different estimations of the exponent a 
according to [24] . 
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Figure 5: Estimations of the exponent a for if 4 theory in two and three dimensions 



9 



(iii) to generalize this mechanism for arbitrary d. 
Surprisingly, this program can be realized and Eq.15 is our main result. Since summation of 
the series gives non-alternating (3(g) (Fig. 4, a), we may conclude that the second possibility 
of the Bogolyubov and Shirkov classification is realized. 



4. "Naive" zero-dimensional limit 

Consider the 0(n)-symmetric (p A theory with an action 

/( n n / " \ ' 

d d x \ E(v^) 2 + \m\ $> a a + §« £ <Pl) 

y a=l a=l \a=l ) 

u = g A € , t = A-d (16) 

in ci-dimensional space; here mo is a bare mass, A is a momentum cut-off, go is a dimen- 
sionless bare charge. It will be essential for us, that the /3-function can be expressed in 
terms of the functional integrals. The general functional integral of (p A theory 

Z ml aM (xi,---,x M ) = / D<p<p ai (xi)<p aa (x2)---tp au (x M )exp(-S{<p}) (17) 



contains M factors of ip in the pre-exponential; this fact is indicated by the subscript M. 

We can take a zero-dimensional limit, considering the system restricted spatially in 
all directions. If its size is sufficiently small, we can neglect the spatial dependence of 
<p(x) and omit the terms with gradients in Eq.17; interpreting the functional integral as a 
multi-dimensional integral on a lattice, we can take the system sufficiently small, so that 
it contains only one lattice site. Consequently, the functional integrals transfer to the 
ordinary integrals: 

Z 2fla M = J dn( P fax ■ ■ ■ <Pa M eX P (~\ m W ~ \ U V~ 

This is the usual understanding of zero- dimensional theory. Such model allows to calculate 
any quantities with zero external momenta. If external momenta are not zero, the model is 
not complete: it does not allow to calculate the momentum dependence. To have a closed 
model, let us accept that there is no momentum dependence at allH. This "naive" model 
is internally consistent but does not correspond to the true zero-dimensional limit of ip A 



5 This point is essential for evaluation of the Z-factor, which is defined in terms of the pair correlator 
G(x — x') — (ip(x)if(x')} in the momentum representation as 

G( ) = 1 = - 

(P> p 2 +ml + E(p, mo) P 2 + m 2 + 0(p i ) ' 

and is determined by the momentum dependence of self-energy. In the described " naive" theory we accept 
Z = 1, since the momentum dependence is absent. 



10 





g 






g* 


t 




Pig) 


t 







Pig) 




a 



Figure 6: a 
of 00/). 



Dependence of g and j3(g) on the parameter £. b — Resulting appearance 



theory. The latter fact is not essential for us, since this model is used only for illustration 
and the proper consideration of the general <i-dimensional case will be given in the next 
section. 

Expressing the /3-function in terms of functional integrals, we obtain it in a form of the 
parametric representation 

n K A K 

9 = 1 - „ , r ^ 2 (19) 



n + 2 Kl 
2n K 4 K " 



n + 2 Kl 



2 + 



K 6 Kq _ -i 

K4.K2 



1 - 



K A K 



(20) 



The right hand sides of these formulas contain the integrals 



K M (t) = V M+n ^d V exp (-t V 2 - ^) , t= (-) 



2 vl/2 



(21) 



obtained from (18) by simple transformations. According to (19,20), the quantities g and 
(3 are functions of the single parameter t; excluding t we obtain the dependence /3(g). 

Investigation of (19, 20) for real t shows that g and (3 as functions of t have a behavior 
shown in Fig. 6, a; combination of these results shows that (3(g) behaves as in Fig. 6, b. We 
see that variation of the parameter t along the real axis determines /3(g) in the finite interval 
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< g < g*, where g* is a fixed point[§ 

g* = — . (22) 
y n + 2 v ; 

To advance into the large g region, we should consider the complex values of t. 

It appears, that in the complex t plane we should be interested in zeroes of the integrals 
i^Af(t). The origin of these zeroes is very simple. There are two saddle points in the integral 
KM(t)i the trivial and nontrivial, 



= , <p c2 = \j-t/2 , (23) 

and K M (t) can be presented as a sum of two saddle point contributions: 

K M {t) = A ie ^ +A 2 e*. (24) 

If these two contributions compensate each other, then the integral can turn to zero. Such 
compensation can be obtained by adjustment of the complex parameter t, and in fact there 
are infinite number of zeroes lying close to lines argt = ±37r/4 and accumulating at infinity 
(Fig. 7). The above saddle-point considerations can be rigorously justified for zeroes lying 
in the large \t\ region. In fact, it is only essential for us that (i) zeroes of K M (t) exist in 
principle, and (ii) zeroes of different integrals lie in different points. 

Now return to the parametric representation (19,20). It appears, that large values of 
g can be achieved only near the root of the integral Ki- If K2 tends to zero, then (19, 20) 
are simplified, 

9 n + 2 Kl ' lA9) n + 2 K$ ' ^ 
and the parametric representation is resolved in the form 

(3(g) = 40, g^oo. (26) 

We see that, indeed, the asymptotic behavior of /3(g) appears to be linear. 



5. General (/-dimensional case 



The same ideas can be applied to the general <i-dimensional case. First of all, the actual 
functional integrals can turn to zero by the same reason. Indeed, the complex values of 
t with large \t\ correspond to complex go with small \qq\ (see Eq.21), and we come to a 
miraculous conclusion: large values of the renormalized charge g corresponds not to large 
values of the bare charge go (as naturally to think0), but to its complex values; more than 

6 Existence of the fixed point g* (obtained previously in [17]) does not mean the existence of a phase 
transition, which is absent for d < 2 due to a finiteness of m 2 . 

7 It is commonly accepted that the bare charge go is the same quantity as the renormalized charge g at 
the length scale A -1 . In fact, these two quantities coincide only on the two-loop level [29] and this relation 
is valid only in the weak coupling region. 
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Figure 7: Zeroes of the integrals Jfjvf(i) in the complex t plane. 

that, it is sufficient to consider the region \g \ <C 1, where the saddle-point approximation 
is applicable. As a result, the zeroes of the functional integrals can be obtained by the 
compensation of the saddle-point contributions of trivial vacuum and of the instanton 
configuration with the minimal action; contributions of higher instantons are inessential 
for \g \ < 1. 

Now we need a representation of the /3-function in terms of functional integrals. The 
Fourier transform of (18) will be denoted as Km after extraction of the ^-function of the 
momentum conservation and a factor I ai ... aM depending on tensor indices: 

(27) 

where J\f is the number of sites on the lattice, and I ai ...a M * s a sum °f terms like S aia2 5 a3a4 . . . 
with all possible pairings. In general, integrals Km{Pi) are taken at zero momenta, and 
only the integral K2 should be known for small momentum 

K 2 (p) = K 2 -K 2 p 2 + ... (28) 

Expressing the /3-function in terms of functional integrals!!, we have a parametric repre- 

8 Definition of the /3-function depends on the specific renormalization scheme. We accept renormalization 
conditions at zero momenta (see Sec. VI. A in [3]). 
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sentation (see [22] for details): 




m K 4 K Q 



K. 



2 



(29) 




, K 4 K 2 (KjKo + g^jgg ~ 2K±K K 2 K 2 
Kl Kl K 2 K 2 -K 2 K 2 



where the prime marks the derivatives over rag. If go and A are fixed, then the right hand 
sides of these equations are functions of only ra , while dependence on the specific choice 
of go and A is absent due to general theorems [3]. 

We see from Eq.29 that large values of g can be obtained near the root of either K 2 , 
or K 2 . If K 2 —> 0, equations (29,30) are simplified, so g and (3 are given by the same 
expression apart from a factor d, 

d/2 K 4 K n , ( K 2 \ d/2 K±Ko 




ki ■ *> = -<{£) -icf- < 31 > 

and the parametric representation is resolved as 

(3(g) = dg, g^oo. (32) 

For K 2 — > 0, the limit g — > oo can be achieved only for d < 4 and we have analogously: 

P(g) = {d-4)g, g^oo. (33) 

The results (32), (33) correspond to different branches of the analytical function (3(g). It 
is easy to understand that the physical branch is the first of them. Indeed, it is well known 
from the phase transitions theory that properties of </3 4 theory change smoothly as a function 
of space dimension, and results for d = 2, 3 can be obtained by an analytic continuation 
from d — 4 — e. According to all available information, the four- dimensional /3-function 
is positive, and thus has a positive asymptotics; by continuity, the positive asymptotics is 
expected for d < 4. The result (32) does obey these demands, while the branch (33) does 
not exist for d = 4 at all. Eq. 32 agrees with the approximate results discussed in Sec. 3 
and with the exact asymptotic result /3(g) = 2g, obtained for the 2D Ising model [9] from 
the duality relation^]. 

6. Strong coupling asymptotics in QED 

The same ideas can be applied to QED. Summation of perturbation series for QED [26J 
gives the non-alternating /3-function (Fig. 8) with the asymptotics (3oog a , where (Fig. 9) 



9 Definition of the /^-function in [9 differs by the sign from the present paper. 
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Figure 8: General appearance of the /3-f unction in QED 

a = 1.0 ±0.1, 000 = 1.0 ±0.3 (34) 

(g = e 2 is the running fine structure constant). Within uncertainty, the obtained /3-function 
satisfies inequality 

0<P{g)<g, (35) 

established in [121 US] from the spectral representations, while the asymptotics (34) corre- 
sponds to the upper bound of (35). Such coincidence does not look incident and indicates 
that the asymptotics 0(g) = g is an exact result. We show below that it is so indeed. 

The general functional integral of QED contains M photonic and 2N fermionic fields 
in the pre-exponential, 

hi,2N = J DADipDip A^{x x ) . . . A llM (x M )ijj(y 1 )i>(z 1 ) . . . ip(y N )i>(z N ) exp ^-S{A,i/j,i/j}^ , 

(36) 

where S{A, ip, ip} is the Euclidean action, 



S{A,ip,ip} = / d 4 x 



- d u A^f + ^(z0 - m + e 4)ij 



(37) 



while eo and mo are the bare charge and mass, and the crossed symbols are convolutions 
of the corresponding quantities with the Dirac matrices. Fourier transforms of the inte- 
grals Im n with excluded ^-functions of the momentum conservation will be referred as 
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K MN (qi,pi) after extraction of the usual factors depending on tensor indices^; and Pi 
are momenta of photons and electrons. 

In general, these functional integrals are taken for zero momenta, but two integrals 
K02 (p) and K2o(q) should be estimated for small momenta: the first is linear in p, and the 
second is quadratic in q, 

K 02 (p) = K 02 + iW , K 20 (q) = K 20 + K 20 q 2 , (38) 

and in fact the tilde denotes their momentum derivatives. 

Expressing the /3-function in terms of functional integrals (see for details), we have 
a parametric representation 

9 = (39) 

-ft- 02 -^20 

2K 02 K( )2 -K> )2 K 02 K 2 2 K 20 \ K 12 K 00 K 02 K 20 j 

where the prime denotes differentiation over tjiq. According to Sees. 4, 5, the strong coupling 
regime for renormalized interaction is related to a zero of a certain functional integral. It 
is clear from (39) that the limit g — » 00 can be realized by two ways: tending to zero either 
Kq2i or K 2 q. For K 02 — > 0, equations (39,40) are simplified, 

K 2 12 K 00 Kf 2 K 00 
9 = - ~ 2 - , P(g) = - ~ 2 ~ , (41) 
^02^20 ^02^20 

and the parametric representation is resolved in the form 

P(g) = g, g^oo. (42) 

For K 20 — > 0, one has 

(3(g) (xg 2 , g^oo. (43) 

Consequently, there are two possibilities for the asymptotics of /3(g), either (42), or (43). 
The second possibility is in conflict with inequality (35), while the first possibility is in 
excellent agreement with results (34) obtained by summation of perturbation series. In our 
opinion, it is sufficient reason to consider Eq.42 as an exact result for the asymptotics of 
the ^-function. It means that the fine structure constant in pure QED behaves as g oc L~ 2 
at small distances L. 



7. Concluding remarks 

10 A specific form of these factors is inessential, since the results are independent on the absolute nor- 
malization of e and m. 
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As should be clear from the preceeding discussion, the conventional renormalization 
procedure defines theory only for < g < g max , where g max is finite. For values g max < 
g < oo, the theory is defined by an analytic continuation, and large values of g correspond 
to complex values of go. Physically, the latter situation looks inadmissible: the S- matrix 
can be expressed through the Dyson T-exponential of the bare action, and Hermiticity of 
the bare Hamiltonian looks crucial for unitarity of theory. 

In fact, a situation is more complicated, as demonstrated by Bogolyubov's axiomatical 
construction of the S- matrix [3]: according to it, the general form of the S-matrix is given 
by the T-exponential of iA, where A is a sum of (i) the bare action, and (ii) a sequence 
of arbitrary "integration constants" which are determined by quasi- local operators. In the 
regularized theory we can set the "integration constants" to be zero, and the 5-matrix is 
determined by the bare action. However, in the course of renormalization these constants 
are taken non-zero, in order to remove divergences. These non-zero "integration constants" 
can be absorbed by the action due to the change of its parameters. As a result, for the 
true continual theory the ^-matrix is determined by the renormalized action, while the 
bare Hamiltonian and the Schrodinger equation are ill-defined. From this point of view 
there is no problem with the complex bare parameters, since the renormalized Lagrangian 
is Hermitian for real g. 

Some problems remain for regularized theory, where the bare and renormalized La- 
grangians are equally admissible and a situation looks controversial. The analogous situa- 
tion was discussed for the exactly solvable Lee model [Hj, which also has the complex bare 
coupling for the sufficiently large renormalized coupling. After the paper [TU] it was gen- 
erally accepted that the Lee model is physically unsatisfactory due to existence of "ghost" 
states (i.e. the states with a negative norm). Quite recently [1] it was found that this point 
of view is incorrect and the Lee model is completely acceptable physical theory. It is a 
key idea of pQ that an analytical continuation of the Hamiltonian parameters to the com- 
plex plane should be assisted by a modification of the inner product for the corresponding 
Hilbert space, 

(/, g) = f f{x)g{x)dx — ► (/, g) G = (/, Gg) , 

and with the proper choice of the operator G the bare Hamiltonian is Hermitian in respect 
to the new inner product (f,g)c- As a result, all states of the Lee model have a positive 
norm and evolution is unitary. The analogous procedure should exist in the present case, in 
order to remove the indicated controversy. In fact, a definition of charge is ambiguous due 
to ambiguity of the renormalization scheme (arising from arbitrariness of " integration 
constants" in Bogolyubov's construction) and complex-valuedness of g has a relative sense 
(see Sec. 5 of [25J). 

The result a = 1 corresponds to one of the really existing branches of the /3-function, 
analytically continued from the weak coupling region. Strictly speaking, we did not prove 
that this branch is physical. This point, together with complex-valuedness of go, casts 
certain doubt on the physical relevance of this result. However, our approximate summation 
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results (Sees. 3, 6), the exact result for the Ising model |9j and inequality (35) for QED give 
the essential evidence that the result a = 1 is physical. 



In conclusion, summation of perturbation series gives the positive /3-function in four- 
dimensional (p 4 theory and QED, while its strong coupling asymptotics is shown to be linear. 
It means that the second possibility in the Bogolyubov and Shirkov classification (Sec. 3) 
is realized, and it is possible to construct the continuous theory with finite interaction at 
large distances. Ej] 
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